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Abstract 

In this paper we consider TV-species nonautonomous competitive systems of par- 
tial differential equations of Kolmogorov type. Under the Neumann boundary 
conditions we give a sufficient condition for the system to be uniformly stable 
and globally attractive. 
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1. Introduction 

In this paper we consider Kolmogorov systems of reaction-diffusion parabolic 
partial differential equations (PDEs) 

du ' 

-g7- = Am, + fi(t,x,ui, . . . ,UN)ui t > 0, x e 0, 1 < i < N. 

From the biological viewpoint this problem represents a model of population 
growth where Ui(t,x) is the density of the i-th species at time t and spatial 
location x € f2, C R n is a bounded habitat, fi(t,x,ui, . . . ,un) is the local 
per capita growth rate of the i-th species. The system is usually considered 
with appropriate boundary conditions, either Dirichlet or Neumann, or Robin 
boundary conditions. In this paper we investigate the system under the Neu- 
mann boundary conditions. We deal with competitive systems. It means that 
the derivatives §^r, 1 < h J ' < N are nonpositive. 

One of the major problems in models of population growth are permanence 
and stability. The authors of many papers focus on finding necessary or sufficient 
conditions for permanence or stability in models which are described by ordinary 
or partial differential equations. Permanence means that any positive solution 
of the model has all its components, for large t, bounded away from zero and the 
lower bound is independent of the solution. Global attractivity states that the 
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difference of any two positive solutions tends to zero as t goes to infinity. In [l| S. 
Ahmad and A.C. Lazer, using the upper and lower averages of a function, found 
sufficient conditions for permanence and global attractivity of the competitive 
N species Lotka-Volterra system of ODEs. In [6j we extended their result for 
the N species competitive Kolmogorov system of ODEs. In terms of the upper 
and lower averages of a function we found sufficient conditions for such a system 
to be permanent and globally attractive. 

In Q the authors considered two species Lotka - Volterra system of PDEs. 
They investigated behavior of this system not only as t goes to infinity, but in 
the "pullback" sence. It means that we are starting with a fixed initial condition 
further and further back in time. They introduced a notion "pullback perma- 
nence" which has an interesting biological interpretation. When two species 
have already been competing for a long time then neither species will have died 
out. Their densities are bounded below in a uniform way, no matter how long 
this population has been running. 

In [2] we investigate the N species competitive Kolmogorov system of PDEs. 
Using the methods of supersolutions and subsolutions for parabolic PDEs we 
found sufficient conditions for permanence in such systems. 

To establish the global attractivity many authors use appropriate Lyapunov 
function in the case of ODE or Lyapunov functional in the case of PDE (see for 
example Q, Q, @|, [§])• In 0, 0| K. Gopalsamy dealt with, respectively, 
periodic and almost periodic solution of a Lotka - Volterra system of ODEs. A. 
Tineo Q considered the nonautonomous Kolmogorov system of ODEs 

uj = / i (t,u)u i l<i<A^, (K) 

where / :Rx — > M. N is a continuous function and is the nonnegative 
orthant ={»£ R N : Ul > 0, 1 < i < N}. He found a sufficient condition for 
the system fKj to be globally attractive. In Q we considered system (|K]) with 
slightly weaker assumptions for the function / than those / in [9( . In [8| Sze-Bi 
Hsu examined Lyapunov functions (or functionals) for various ecological models 
which take form of ODE systems (or reaction-diffusion PDE systems). First 
the author constructed Lyapunov functions for the predator prey system and 
then constructed Lyapunov functionals for the corresponding reaction-diffusion 
PDE system. Sze-Bi Hsu considered the models in which the reaction term is 
independent on t and x. 

In A. Tineo and J. Riviero dealt with the two species nonautonomous 
competitive Lotka-Volterra system with diffusion. They resigned from the Lya- 
punov function (functional) to establish of attractivity of this system. They 
used a powerful tool which is the iteration monotone methods. 

Usually the same conditions which guarantee permanence in the ODE model 
of population growth ensure that the model is globally attractive ( see for ex- 
ample 0). But when we try to prove that a system of Kolmogorov type of 
PDEs is globally attractive then some difficulties occur. Firstly, we can prove 
that the system is globally attractive only under the Neumann boundary con- 
ditions. Secondly, besides the conditions which guarantee the permanence we 



2 



need additional inequalities which contain the parameters which occur in the 
definition of permanence. Even in the case of two species problem it is a diffi- 
cult to obtain sufficient conditions which guarantee both permanence and global 
attractivity of Kolmogorov system of PDEs. One can try to introduce a partial 
ordering. The idea is that we can compare solutions of our system of PDE with 
the functions which are independent on x in a way that when we apply partial 
ordering our solutions of PDE are between the functions independent on x. A 
difficulty occurs when we let the parameter which is a spatial location to be 
arbitrary. 

In this paper we find sufficient conditions for N - species system of PDEs 
of Kolmogorov type under the Neumann boundary conditions to be globally 
attractive. Conditions which we obtained are stronger than those guaranteeing 
perrmanence. Assumptions for a function / in this paper is the same as in 

This paper is organized as follows. 

In Section 2 we present basic assumptions and definitions. 

In Section 3 we formulate the main theorem of this paper. By Lyapunov 
functional some sufficient conditions are obtained to guarantee attractivity and 
uniform stability Kolmogorov system of PDEs. 

In Section 4 we give some remarks. 



2. Preliminaries 

We consider a nonautonomous system of PDEs of Kolmogorov type 
du ' 

— -± = Aui + fi(t,x, ujl, . . . ,u N )ui, t>0, xeil, i = l,...,N 
9t (R) 

Bu t = 0, t > 0, x € dQ, i = 1, . . . ,N, 

where f2 C K n is a bounded domain with a sufficiently smooth boundary <9f2, A 
is the Laplace operator on Q, 

dx\ dx 2 n ' 

and B is the boundary operator of the Neumann type 

Bui = on 90, 

ov 

where v denotes the unit normal vector on <9il pointing out of f2. 
Now we make the following assumptions 

(AI) (see [2]) /»: [0, oo) x fi x [0,00)^ -» R (1 < i < N), as well as their 
first derivatives dfi/dt (1 < i < N), dfi/duj (1 < i,j < N), and dfi/dxk 
(1 < i < N, 1 < k < n), are continuous. 

(AI) guarantees that for any regular initial function Uq(x) — (ttoi, • • • > u on), 
x € O where uoi > (1 < i < N, x £ fi) there exists a unique maximally 
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defined solution u(t,x) — (ui(t, x), . . . , Uff(t, x)) of fR|, (t,x) G [0,r max ) x Cl, 
where r max > 0, satisfying the initial condition u(0,x) — Uq(x). Moreover, the 
solution of system |R| is classical. 

(A2)(see [2]) Functions [[0,oo) x fi 3 (t,x) i-> f t (t,x,0,. . . ,0) G R], 1 < i < 
JV, are bounded. 

Definition 1. j4 solution u(t,x) — (ui(t, x), . . . , Ujj{t, a;)) of system (|R|) is 
positive if Ui(t, x) > /or all i = l,...,N, x G O, t > 0. 

Definition 2. System (|R| is permanent i/ £/iere ea;isi positive constants 5, 5 
such that for any positive solution u(t, x) = (u±(t, x), . . . , uat(£, a;)) £/iere exists 
T = T(u) > i/ie property 

5 < Uiit, x) < 5 
for allt>T, x i = l,...,N. 

Definition 3. System (|R| is globally attractive i/ any too positive solutions 
u(t, x) = (ui(t, x), . . . , U]y(t, x)) and v(t, x) = (vi(t, x), . . . , Vpf(t, x)) of system 
(|R|) satisfy 

lim x) — a;)) = 
for 1 < i < A, uniformly in x G fi. 

Definition 4. System is uniformly stable on [ii,oo), ti > if for ev- 
ery e > i/iere is a (5 > so </ia< /or any too positive solutions u(t, x) = 
(tti(t, a;), . . . , ttjv(t, a;)) and v(t,x) = (i>i(£, a;), . . . , «j\r(t, a:)) i/ £2 G [ii,oo) and 
|n(t 2 , 2;) — u(<2) aOH < ^ i/ien \\u(t, x) — v(t, x)\\ < e for t > t 2 - 



The next assumption is 

(A3)(see 0|) There exist b u > such that Qi±(t,x,u) < -b u for all t > 0, 
x en, ue [0,oo) N , 1 <i < N. 

Define 



B(e) := 




x ■ 


■ X 


0,- +c 


. hi 






. 2aw 



(A4)(see 0) dfi/duj, 1 < i,j < N, are bounded and uniformly continuous 
on each set [0, 00) x Cl x B. 

For 1 < i,j < N define 

( dfi 

bij(e) :— sup< — — — -(t,x,u) : t > 0, x G Q, u G J5(e) 
and 6ij(0) := bij. 
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Assumptions (A3) and (A4) imply that 6y(e) > 0, 1 < i, j, < N, and bu(e) > 
0, 1 < i < N, whereas it follows from (A5) that by (e) < oo, and lim c ^ + hj (e) = 
kj, for 1 < i,j < N. 

Recall that in |2| we found average conditions for system |R| to be per- 
manent. The average conditions for system fR| under the Neumann boundary 
conditions are 

mlft]> ±hfm, (AC) 

where 

1 f* 

m[fi] := liminf / min/^r, x, 0, . . . , 0)dr, 

t-s-Hx> t — s J s xen 

and 

1 f* 

M[fi\ := hmsup / max/j(r, x, 0, . . . , 0)dr. 

t— s— >oo i S J s i£f! 



3. Main Theorem 

In this section we state the main theorem of this paper and with the help of 
a Lyapunov functional we prove attractivity and stability of system ([R"j) . 

Theorem 3.1. Assume (Al) - (A5) and (AC). Let 5, 6 > be such that for 
each 1 < i < N 

N 

3=1 

Suppose that u(t, x) — (u\(t, x), . . . , ujv(£, x)) and v(t, x) ~ (vi(t, x), . . . , ujv(i, x)) 
be two positive solutions of system (|R|) sitc/i that 5 < Ui(t,x),Vi(t,x) < d for 
sufficiently large t and all x £ O. Then there exist positive constants Z , 7 and 
T > smc/i </ia< /or i = 1, . . . , N and t >T one has 



N N 

sup|w 4 (i,x) - ^(i,x)| < sup|ui(r, a;) - Uj(T,a;)| • exp (— 7(t - T)), 

i=i ^en i=1 xen 



In particular system (|R)) is globally attractive and uniformly stable. 

Proof. Fix two positive solutions u(£, a;) = (ui(i, a;), . . . , u^(t, x)) and x) = 
(v\(t,x), . . . ,vpf(t, x)) of system (|R]). By assumption (A5) we can take e > 
such that 

N 

<Li hi > ^ Sj hj (e) for all i = 1, . . . , N. 

3=1 



Let to > be such that u(t, x), v(t, x) G B(e) and 5 < Ui(t,x),Vi(t,x) < S for 
all t > to, x G O, 1 < i < N. By a result in matrix theory (see [7|, Section 5]) 
we know that there exist a%, . . . , a at > such that 



JV 



a i<Lku > ^ a j$bj l (s), 1 < i < N. 



u 2 (t, x) 



Vi(t,x) 



Denote 

N 

9(i) = '^^a i O i (t) 1 where Qi(t) :— sup In 

i=i x ^ 

Fix 1 < i < N. We prove that there exists T > to such that 

JV 

D+Qi(t) < -Sb^i^ + Sj^hA^M t > f. 

i=i 



(3.1) 



(3.2) 



(3.3) 



where D + Q denotes the upper derivative of O. 

We take for T any number > i . For t >T denote Cl^(t) to be the set of the 
points a;W g such that the supremum of the function 



ThenCi^(t) = ^ ) (i)u^ ) (t), n^(t)nfiW(t) 



>(*)/ 



ln Hi(*4 



is realized. 



where 



f)W(t) = {a;W g nW(t) : m(t,a;«) > «i(t,a;W)} 



and 



fiW(i) = {a; w G flW(t) : U! (t,i w ) < «i(t,a;W)} 
Therefore by (JOJ) 



-D + Oi(i) < max< sup 



(3.4) 



Denote m := &i(t). Now we consider four (not mutually exclusive) cases. 

Case one. G finA^(t). 

By the definition of n it follows that 

Ui(t,x) < e^Vifax) for all x G and U;(t, x w ) = e^v^t, x^), (3.5) 

and 



d_ 

dt 



In- 



Ui(t,xW)7 m(t,a;W) 



(t,iW) 



iti(i,xW) Vi[t,xW) 
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By (EU), 
Hence 

dt \ Vi(t,xW)J 

Case two. X® G fin#(t). 

In a similar manner we prove that 

Jt ( ln ^xWj) ~ Mt > X(i) ' V(t > X(t))) ~ Mt ' Xit) ' U{t > X(l))) - 
Case three, x« £ dfl D n^(t) 

This implies that Ui(t,x) < e**Ui(i, a^) for all x € f2 and Uj(f,a;W) = 
e^v^t, x^). Denote 5(2;) := Ui{t,x) — e Xi v%{t, x). The restriction g\en of g to 
<9f2 has largest value at the point x^ l \ consequently V(g\da)(x^ 1 ') = 0. By the 
definition of g it follows that §2 = |H - = 0. Hence Vff(a?W) = 0. 

As (7 attains its maximum at a;W, this implies that Ag(x^') < 0. Therefore 

d f ln Ui{t,xW)\ _ x«) £vi(t,xW) 



dt V Vi(t,xW)J u»(i,a;W) 

< /((*, x« , u(t, - /<(*, x« w (t, iW)). 

Case four, <E OU#(i) 

In a similar manner we prove that 

^(^ ^S ) </i(*,a: W ,t;(t, a; W))-/ i (t, a; W,u(t, a; W)). 



The cases 1-4 with (13.41) give 

i=l 

sup{/ i (i J i( < ^u(t,a:W))-/ < (t,a:( < \t>(t,s«)) en^}}. 
By assumptions (A3), (A4) and (A5), it follows that 

f^x^ , Ul (t,x^), . . . ,u N (t,x^)) ~ f^v^x^)), . . . ,v N (t,x^)) 

N 

< -b H ( Ui (t, x« ) - (t, + £ &« (e) I % (*, 1 W ) ~ ^ (W°)| 



7 



for each x« G 

By assumptions (A3), (A4) and (A5), it follows that 

fi{t, x^), Ul (t, x«), . . . , Ujv (t, X W) _ f t ( t , Vl (t, iW)), . . . , ujv (i, x«)) 



< 



(«i (t , ) - u< (t, x« ) ) + £ Ey (e) | ^ (t, x« ) - Uj (t, x« ) | 



3=1 



for each x« G 0^(4). 
Hence 



A' 



D + 9j(t) < max{ sup {-^M*, x®) ~ «,(*, x«)) + £ E«(e)l«j(*> x«) - Wj (t, x«)|}, 



N 

sup («< (t , x« ) - Ul (t , i W )) + ^ by (e) | Wi (i , i W ) - u,- (t, x W ) | } } 

arWefiW 3=1 

3#« 

Since 5j < x), i>j(i, x) < i?i for alH = 1, . . . , N and for sufficiently large 
t, using the mean value theorem we have that 



= \uj(t,x) - Vj(t,x)\ < 



In 



Uj(t, x) 



Vj(t,x) 



< -\uj(t,x) -Vj(t,x)\ (3.6) 



for all t > T, x G O, 1 < j < N. 
Therefore 



L> + e t (i) < max{sup{-£6, 
sup{-Sb it 

Hence we proved p.3[) . 



In 



t (t,x^) 



N 



In 



«i(t,a:W) 



3=1 

N 



ln %(M (,> ) 



sup 

£Cgf2 



In 



Uj(t, x) 



Vi{t,x) 



3=1 
3#» 

N 



hi 



Uj(t,xM) 



X« G 



«,•(*, x«) 



: xW G Q« 



>} 



3=1 



In 



%(t, x) 



t>j(*,x) 



f > T. 
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By l|3.ip . it follows that 

N / N 



i=l ^ 3=1 y 

N N N 

= -JT i a i 5b u Q i {t)+-SY i ^M^M 

i—1 j — 1 i— 1 

N N 

»=i j=i 

3*7^ 

AT 

<-e^ei(t), t>T, 



where 



f W 

-~ ^ i'=l 



3 = 
3¥i 



Hence by ([33j and ([3J|) . 

AT AT 



^ + e(t) ^^aiU+GiCt) < -e5^0i(i) < -6(t), t>f (3.8) 



(3.7) 



where a* = maxja^ : 1 < i < N}. Therefore 

e(t)<0(f)exp(-^(t-f)) (3.9) 



for all i > T. 

By it follows that 

AT N 

e(t)>a*5^e i (t)>^5^sup|« i (t,a:)-« i (i,x)|, t>T, (3.10) 

and 

6(f) < a*^6i(t) < —J2sup\ Ui (f,x) - Vi (f,x)\, t>f (3.11) 
i=i - i=i*efi 

where a* = min{aj : 1 < i < TV}. 

From (jnU), (|3~TU| and (gUP we have that 

TV N 

sup|ti»(i,a:)— x)| < sup|ui(T, x)-v l (f, x)|-exp (-j{t - T)), i > T, 

i=i i=i ^en 

(3.12) 
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where Z — j±- and 7 = -% [Z and 7 we can take independently on the solutions 
u(t, x) and vit, x)). 

By (|3.12| it follows that the system (jRj is globally attractive and uniformly 
stable. 

□ 



4. Remarks 



Remark 1. Note that conditions (|3.ip contain <5 and (5 which appear in 
definition of permanence of system (jR|. In [2| we showed that if we change 
conditions (jACj) to slightly stronger that is 

N T _ 

-[/*]> E^p ( AC ') 
3=1 

then we can estimate 8 and (5 in terms of parameters of system (|Rj) . By dissipa- 
tivity of system JR| it follows that we can take maxi^^jvjf^} as S. Concerning 
the lower bound (see [2|, Theorem 4.1]) as 5 we can take 

1 N T - 



min {= L( a ._y« } 



3=1 ~" 
when 

N T - 

3=1 -v 

for aU 1 < i < N. 

If we define S_ and 5 as in the above condition (12.11) , take a form 



m > E 



x ' b M M[f j \ bjjM[fi\ 
en — h aj 

hi hi 



Remark 2. 



Definition 5. System is uniformly asymptotically stable on [ii,oo) if it 
is uniformly stable and if there is a r > so that for every e > there is a 
T(e) > so that for any two positive solutions u(t,x) = (ui(t, x), . . . , u/v(i, x)) 
and v(t, x) = {v\{t, x), . . . , Ujv(i ; x)) if ta > ii a«d 1 1 u 2 5 x) — vfo, x) \\ < r then 
\\u(t,x) -v(t,x)\\ <e fort> t 2 +T(e). 

Note that from attractivity and uniform stability of system (|R]) it follows 
that system (JRj) is uniformly asymptotically stable. 
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